Superflow in a toroidal Bose-Einstein condensate: 
an atom circuit with a tunable weak link 
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We have created a long-lived (~ 40 s) persistent current in a toroidal Bose-Einstein condensate 
held in an all-optical trap. A repulsive optical barrier across one side of the torus creates a tunable 
weak link in the condensate circuit, which can affect the current around the loop. Superflow stops 
abruptly at a barrier strength such that the local flow velocity at the barrier exceeds a critical 
velocity. The measured critical velocity is consistent with dissipation due to the creation of vortex- 
antivortex pairs. This system is the first realization of an elementary closed- loop atom circuit. 
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Quantum fluids can exhibit properties such as long- 
range coherence and superfluidity that make them useful 
for constructing sensors and other devices. For example, 
superconducting quantum interference devices (SQUIDs) 
are sensitive magnetic field detectors pQ, and superfluid 
He circuits have been used to detect rotation |2l[3]. Ultra- 
cold atomic-gas analogs of electronic devices and circuits, 
or "atomtronics" have been proposed including diodes 
and transistors [4]. Of particular interest is the realiza- 
tion of an atomic-gas SQUID analog. SQUID circuits 
have been realized with either tunnel or weak link junc- 
tions [H 111 |6]. In atomic Bose-Einstein condensates, 
Josephson junctions have been demonstrated only be- 
tween adjacent wells [71 [8]. Here we present the first 
implementation of a non-trivial, closed- loop atom circuit, 
and show that it is possible to control the current at the 
single-quantum level by changing the strength of a weak 
link. This is an essential step toward realizing an atomic 
SQUID analog. 

Superfluids flow without dissipation if the flow veloc- 
ity is below a threshold determined by the lowest en- 
ergy excitations allowed for the system [9]. In a ho- 
mogeneous condensate the lowest energy excitations are 
phonons [10], and the Landau critical velocity is the 
speed of sound [11] . Real systems are finite and therefore 
inhomogeneous; consequently the lowest energy excita- 
tions can be vortex-like [12^, and dissipation can occur 
at velocities well below the sound speed [13] . Dissipation 
involving vortex-like excitations has been previously ob- 
served in experiments with liquid He [HI [15] , supercon- 
ductors [16] , and in a simply-connected condensate [17] . 

The critical velocity in simply-connected condensates 
has been measured previously by moving a defect created 
by a localized optical potential [TTHCT. When the veloc- 
ity of the defect was high enough, excitations and heating 
were observed. In contrast to this earlier work, we create 
a quantized, persistent flow around a multiply-connected 
(toroidal) condensate, and study the decay of that flow 
in the presence of a stationary barrier, as a function of 
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FIG. 1: Experimental configuration, (a) Schematic of the 
toroidal optical dipole trap formed at the intersection of two 
red-detuned beams: a horizontal "sheet" beam, and a vertical 
Laguerre- Gaussian beam (LGj) with a ring-shaped intensity 
maximum. A pulsed pair of Raman beams (large downward 
arrows) co-propagating with the LG trapping beam creates 
circulation in the condensate, (b) Energy level diagram for 
the Raman transition: \F — 1, mp — —1) -^\F —1^ mp — 0). 
One Raman beam carries h orbital angular momentum per 
photon (LGJ), the other carries none (Gaussian); the con- 
densate is transferred to a quantized (/=1) circulating state, 
(c) False-color absorption image showing the normalized col- 
umn density of a condensate in the trap, viewed from above. 
Arrows: Raman beam polarizations and magnetic bias. 



barrier height and condensate atom number. 

In previous experiments [20 , we created persistent cur- 
rents in a harmonic magnetic potential pierced by a re- 
pulsive optical potential. Relative drift between these 
potentials limited the flow lifetime to ~ 10 s. This moti- 
vated the construction of an all-optical trap which sup- 
ports persistent currents for up to 40 s, and allows us to 
carefully study the stability of superflow. 

To create a toroidal condensate, 3^5'i/2 |F = 1, m^? = 
— 1) ^^Na atoms are cooled almost to degeneracy in a 
magnetic trap and then transferred into an optical dipole 
trap created by the intersection of red-detuned (1030 nm) 
"sheet" and "ring" beams (Fig.jl^). The horizontal sheet 
beam has a vertical (horizontal) 1/e^ half- width of ~ 
9 /im 400 /im), and provides vertical confinement. 
The vertical ring beam is Laguerre- Gaussian (LGj), and 
confines the condensate to its ~ 20 jiia radial intensity 
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peak, generating a toroidal potential minimum. With 
the atoms in the optical trap, the beam intensities are 
ramped down to force evaporative cooling. At the end of 
the ramp, the trap depth is ~ 700 nK, with trap frequen- 
cies uOz/"^^ = 550 Hz (vertical) and ujrl'^^ = HO Hz (ra- 
dial). This produces a toroidal condensate of up to 3x 10^ 
atoms with a chemical potential /io of up to /z • 1200 Hz, 
and temperature < 10 nK (no discernible non-condensed 
fraction). The azimuthal variation of the potential min- 
imum is less than h • 100 Hz, as shown by the smooth 
condensate density profile in Fig. [T]^. 

The condensate is initially nonrotating [21] . Superfluid 
circulation around any closed path must be quantized, 
such that the wave function has a 2i\l phase winding 
{I gZ). We create circulation by transferring quantized 
angular momentum from optical fields during a Raman 
process [22]. The co-propagating Raman beams, detuned 
2.3 GHz below the D2 transitions, are in two-photon reso- 
nance with the |1, —1) ^ |1, 0) transition (Fig. [l]3). They 
have orthogonal linear polarizations, parallel and perpen- 
dicular to the horizontal magnetic bias field (Fig. 
The nonlinear Zeeman shift from the 0.5 mT field ap- 
plied during the interaction prevents coupling to |1, 1). 

The angular momentum change of the condensate is 
determined by the spatial mode of the Raman beams. 
With one beam Gaussian, and the other in an LGj spatial 
mode carrying h orbital angular momentum per photon, 
the Raman process coherently transfers the condensate 
to the 1=1 circulating state [2ai[22l[23]. With good mode 
matching and an optimized Raman 7r-pulse (^ 100 /is), 
we achieve a minimum transfer efficiency of 90%, with 
only a few percent atom loss due to spontaneous scatter- 
ing. Residual atoms in |1, — 1) after the Raman pulse are 
quickly removed from the trap by transferring them to 
1 2, —2) with a microwave pulse, then ejecting them from 
the trap with resonant imaging light (see below). 

Circulation is detected by releasing the condensate 
from the trap and imaging the density distribution af- 
ter several milliseconds time-of-ffight (TOF) [24]. If the 
condensate is not rotating, the central hole closes after 
a short time. When a rotating condensate is released, 
the angular velocity of the ffow prevents complete clo- 
sure. The persistence of a central hole after sufficiently 
long TOF is the signature of circulation in the ring (see 
Fig[2|b) insets). The apparent size of the hole at a given 
time after release is related to the the azimuthal flow 
velocity prior to release and the velocity of the mean- 
field-driven inward expansion. For a rotating condensate 
released directly from our narrow annular trap, the hole 
size is below our imaging resolution for experimentally 
accessible TOFs (< 15 ms). To make the signature of cir- 
culation visible earlier, we first adiabatically reduce the 
ring beam intensity by 90% over 100 ms, then release the 
condensate suddenly (< 1 /is). We use this procedure, 
followed by 6 ms TOF, to detect circulation. 

The Raman beams used to create circulation also cause 
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FIG. 2: (a) Flow survival as a function of chemical potential, 
/io, for two barrier heights: = 650 Hz (upper, blue), and 
P/h = 780 Hz (lower, red). Presence or absence of flow for 
a single condensate is shown by closed circles. Open circles 
are the average of data within the bins (vertical lines), rep- 
resenting the flow survival probability (Paow) of each bin. A 
critical chemical potential fic for stable flow is found from a 
sigmoidal fit (solid lines) to the data for each /3. Inset: In situ 
absorption image of a condensate near /ic {1^0/ h = 870 Hz, 
P/h = 650 Hz), (b) Values of /ic at different /3, determined 
by flts as described in (a) . The open circles correspond to the 
data in (a). The vertical error bars reflect the width of the 
sigmoidal flt, ±2/i^^. The horizontal error bars are the 20 Hz 
uncertainty in calibrating /3. The dotted line is a linear flt to 
the data, with slope 1.6(2). Insets: typical TOF absorption 
images showing the presence (top left) and absence (bottom 
right) of circulation. 



small-amplitude oscillations in the radial density proffle, 
due to small dipole forces and atom loss. These oscilla- 
tions have no observable impact on the stability of the 
circulation, and damp out after < 0.5 s. We add a wait 
time > 3 s after the Raman transfer to ensure complete 
damping. The circulation is extremely robust, and con- 
tinues until losses due to background collisions reduce /io 
to the level of the nonuniformities in the trap. For a 30 s 
vacuum-limited 1/e condensate lifetime, /io remains high 
enough for ffow to survive up to 40 s. 

After the > 3 s wait time, we insert a barrier into the 
path of the ffow and study the stability of the circulation. 
The repulsive barrier is created with a blue-detuned (532 
nm) laser beam focused to an elliptical spot. The major 
axis (15 /im 1/e^ radius) is aligned in the radial direction 
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of the toroid, and the minor axis (4.3 fim radius) 
is parahel to the flow, and exceeds the bulk condensate 
healing length = h/^/2rn^I < 1 /im). The barrier de- 
pletes the local density of the condensate, n, as seen in 
the inset in Fig. [2|a). The reduction in density increases 
the local flow velocity (roughly v ocl/n). Lowering the 
density also lowers the local interaction energy, /i/ ex n, 
decreasing the local sound speed. To study flow stability, 
we ramp up the barrier intensity over 100 ms to a cho- 
sen barrier height /3, hold for 2 s, then ramp back down 
in another 100 ms. The presence or absence of flow is 
then detected in TOF as described above. This proce- 
dure is repeated many times for the same /3, varying the 
total number of atoms (by varying the initial condensate 
number and/or the wait time) until the range of atom 
number is well-sampled. We then change /3 and repeat 
the procedure. If the barrier is not applied, the flow al- 
ways survives, so we can attribute the decay of the flow to 
the effect of the barrier. Separate measurements indicate 
that the flow decays in < 100 ms. 

The analysis of flow stability depends on in situ ob- 
servations of the condensate density profile in the pres- 
ence of the barrier, and from TOF images after the bar- 
rier has been removed. From TOF images we determine 
whether the flow survived [insets Fig. |2|b)], and mea- 
sure the condensate atom number, N. For an annular 
condensate with a Thomas- Fermi profile, the chemical 
potential /io = hUj^Jnj^ • {Nas/R) where uj = ^JUj^uj^^ 
as is the s-wave scattering length, and R is the radius 
of the ring. This calculation does not include small cor- 
rections (~ 6%) due to the azimuthal nonuniformity of 
the potential minimum and displacement of atoms from 
the barrier region, corrections which are less than the 
systematic uncertainty in determining /io (^ 10%). 

We calibrate /3 by taking in situ images of the conden- 
sate and measuring the reduction in column density at 
the location of the barrier (see Fig.[2|i inset). Due to the 
high optical depth (up to 10), we use a partial transfer 
imaging technique [25l [26], in which a precise fraction 
(ranging from 15-40%) of the atoms is transferred to the 
|2, — 1) state using a microwave pulse, then resonantly im- 
aged on the S1/2F = 2 P^i^F = Z transition. The lo- 
cal interaction energy jii can be found from the measured 
column density n. For data where jii < hwzi we assume 
the axial density profile is that of the harmonic oscillator 
ground state, with jii = [Sir • {hujz){h'^alh'^ /m)]^^'^ . For 
data where /i^ > huz^ we assume a Thomas- Fermi profile, 
with /ii = [97r^/2 • {huJz)'^{^'^CLln'^ /m)]^^^ . When measur- 
ing the column density at the barrier, we correct for loss 
of contrast due to the imaging resolution, which reduces 
the apparent depth of the density depletion by ~ 15%. 
We take f3 to be fiQ — fii. 

The solid circles in Fig. |2|a) show the flow survival or 
decay for single experimental runs, plotted against /io foi" 
that run. The open circles are the average of the solid 
circles within the bins shown. The upper plot (blue) is 
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FIG. 3: Critical flow velocity, Vc, above which the circulation 
becomes unstable, for each of the barrier heights in Fig[2|b). 
The value of Vc is shown normalized to the effective sound 
speed at the barrier (ceff), and is plotted as a function of Ceff . 
The horizontal error bars are the estimated uncertainty in 
Ceff. The vertical error bars are the combined experimental 
uncertainty in Vc and Ceff. The measured ratio Vc/cefi ~ 0.6, 
and is independent of Ceff to within the experimental uncer- 
tainty. Vortex-like excitations are expected to occur in this 
system at and above a velocity vf, where vf < Ceff [13]. The 
gray band indicates estimated upper/lower bounds (see text) 
on VF/cef£, using Feynman's approximate expression for vf- 



for P/h = 650 Hz (upper); the lower (red) shows f3/h = 
780 Hz. At low /io, the flow is arrested by the barrier. 
At high /io the flow survival probability becomes unity. 
In between, the survival probability increases from zero 
to one over a narrow critical region. We characterize this 
critical region for each (3 by fitting a sigmoidal function 



P(/io) = 1/(1 + e 



(/ic-/io)//i^ 



) to each unbinned data set. 



where the parameters /ic and fi^ are the critical chemical 
potential and the sigmoidal width respectively [27 . The 
observed width is consistent with observed shot-to-shot 
variations in the trapping potential. 

Figure [2|b) shows the values of /ic extracted from the 
fits of the data for seven different p. Over this range, 
/ic increases approximately linearly with /3, with a slope 
greater than unity. The functional dependence and slope 
are determined in a non-trivial way by trap geometry and 
the condition of quantized circulation around the ring. 
The experimental results are consistent with expectations 
for our geometry. 

The physics behind Fig. |2] is more apparent when the 
data is recast in terms of flow velocity and sound speed 
at the barrier. The barrier thickness is greater than ^, so 
we expect the flow to become unstable when the velocity 
in the barrier region exceeds some local critical velocity 
Vc. The flow velocity at the barrier cannot be determined 
just from /io and /3. The requirements of quantized cir- 
culation (global), and flow conservation (local), make it 
necessary to self-consistently calculate the velocity distri- 
bution around the entire ring. We do this by integrating 
the in situ column density radially to make a ID approxi- 
mation of the density profile, then solving for the velocity 
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distribution of an / = 1 circulation state. 

The critical velocity is determined by the lowest en- 
ergy excitations allowed for the system [9]. For phonon- 
like excitations in the ring, that velocity should be ap- 
proximately the local sound speed in the barrier re- 
gion ^ . We make an initial estimate for the critical 
velocity from the local interaction energy at the peak 
of the barrier, q = ^J [iijm. However, the inhomoge- 
neous (nearly parabolic) radial density profile lowers the 
effective sound speed to Cgfr = q/ a/2 for waves traveling 
azimuthally along the annulus [29] . Figure [s] shows the 
observed critical velocity normalized to Ceff , as a function 
of Ceff. As seen in previous work with finite inhomoge- 
neous atomic condensates p!7tfl9l . the observed critical 
velocity is less than the sound speed. For all tested val- 
ues of /3, 'L'c/ceff ^ 0.6 and is independent of Ceff to within 
the experimental uncertainty. 

In this experiment, flow is confined to a narrow, flat- 
tened channel, raising the possibility that vortex-like ex- 
citations are responsible for the observed critical veloc- 
ity. Numerical simulations [30] with a model conden- 
sate similar to ours, but in an / = 8 circulation state, 
showed vortices traversing the barrier region when the 
barrier was raised above a critical level. This suggests 
that for our / = 1 circulation state, a similar decay mech- 
anism could be at work. For vortex-like excitations in 
our quasi-2D geometry, the (Feynman) critical velocity 
vp can be estimated from energetic arguments [13 to be 
vf = (h/md) ln((i/a), where d is the channel width, and a 
is the vortex core size. We take d to be the Thomas- Fermi 
width, and a the healing length, both calculated for the 
barrier region. Both d and a depend on Ceff via the inter- 
action energy /i^ . The grey band in Fig. |3] is an estimate 
of the probable value of vf/cq{^ with q > Cgff > q/>/2. 
While this calculation is in surprisingly good agreement 
with our data, a more complete model including geomet- 
ric factors is needed to accurately calculate the energy of 
a vortex-antivortex pair in the barrier region. 

We have presented the first realization of a closed 
atomtronic circuit, demonstrating precise control both in 
inducing and arresting superfluid flow. We have clearly 
identified the critical velocity where flow stops, and our 
observations are in agreement with theoretical predic- 
tions in which vortex-antivortex excitations are the decay 
mechanism for the system. In future work, we plan to 
investigate the role of barrier geometry, condensate tem- 
perature, and dimensionality in determining the critical 
velocity and decay mode. In addition, rotating a bar- 
rier around the ring (oscillating it azimuthally) would be 
analogous to magnetically biasing (driving an AC current 
in) a SQUID. The present work constitutes a significant 
step toward realizing such an atomic SQUID analog. 
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